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We investigate the influence of a temperature-dependent shear viscosity over entropy density ratio 
r;/s on the transverse momentum spectra and elliptic flow of hadrons in ultrarelativistic heavy-ion 
collisions. We find that the elliptic flow in -y/ijvjv = 200 GeV Au-|-Au collisions at RHIC is dominated 
by the viscosity in the hadronic phase and in the phase transition region, but largely insensitive 
to the viscosity of the quark-gluon plasma (QGP). At the highest LHC energy, the elliptic flow 
becomes sensitive to the QGP viscosity and insensitive to the hadronic viscosity. 

PACS numbers: 25.75.-q, 25.75.Ld, 12.38.Mh, 24.10.Nz 



Ultrarelativistic heavy-ion collisions at the Relativistic 
Heavy-Ion Collider (RHIC) and the Large Hadron Col- 
lider (LHC) produce a hot and dense system of strongly 
interacting matter The subsequent expansion of the 
created matter has been shown to exhibit a strong de- 
gree of collectivity which reveals itself in the transverse 
momentum {px) spectra of finally observed hadrons. In 
particular, the observed large azimuthal anisotropy of the 
spectra, quantified by the so-called elliptic flow coefficient 
V2, has been interpreted as a signal for the formation of a 
quark-gluon plasma (QGP) with very small viscosity in 
heavy- ion collisions at RHIC Q . 

A first indication for the small viscosity of the QGP 
was the agreement between RHIC data and hydrodynam- 
ical calculations in the perfect-fluid limit, i.e., with zero 
viscosity Q. An analysis of the elliptic flow at RHIC in 
the framework of relativistic dissipative hydrodynamics 
was performed in Refs. These works indeed indi- 

cate that the shear viscosity to entropy density ratio, rj/s, 
has to be small in order to keep the agreement between 
the hydrodynamic simulations and experimental data. 

Presently, most hydrodynamical simulations assume 
a constant, i.e^ temperature-independent ry/s. It has 
been claimed [q that, in order to describe elliptic flow 
data, this value cannot be larger than 2.5 times the lower 
bound 77/s — 1/Att conjectured in the framework of the 
AdS/CFT correspondence 0|- A constant ij/s is, how- 
ever, in sharp contrast to the behavior observed in com- 
mon liquids and gases, where rj/s has a strong temper- 
ature dependence and, typically, a minimum near phase 
transitions. A similar behavior of rj/s is expected for 
finite-temperature matter described by quantum chro- 
modynamics (QCD) near the transition from hadronic 
matter to the QGP (the QCD phase transition) Q. 

A natural question then is whether the temperature 
dependence of 77/s has an effect on the collective flow of 
hadrons in heavy-ion collisions. In this work, we inves- 
tigate this question in the framework of relativistic hy- 
drodynamics. Wc consider a temperature-dependent 77/s 



with a minimum near the QCD phase transition, and 
compare the results with those obtained for a constant 
77/s in either the hadronic phase, or the QGP phase, or 
both phases. Note that we do not attempt a detailed fit 
to the data in order to extract 77/s. Rather, we are inter- 
ested in the qualitative effects of different parametriza- 
tions for 77/s on hadron spectra and elliptic flow. 

Concerning the elliptic flow in Au-|-Au collisions at 
RHIC, we find little difference whether rj/s is constant in 
the QGP phase or strongly increasing with temperature. 
In contrast, the elliptic flow values are highly sensitive 
to whether we use a constant or temperature-dependent 
77/s in the hadronic phase, corroborating the findings of 
Refs. On the other hand, we find that the sen- 

sitivity of the elliptic flow to the values of 77/s in the 
high-temperature QGP increases with increasing collision 
energy, while the sensitivity to the hadronic viscosity de- 
creases. At the highest LHC energy, the above conclusion 
for RHIC energies is reversed: the finally observable el- 
liptic flow is dominated by the viscosity of the QGP and 
largely insensitive to that of the hadronic phase. 

Fluid dynamics is determined by the conservation of 
energy, momentum, and charges like baryon number. 
Here, we are interested in the collective flow at midra- 
pidity in heavy-ion collisions at RHIC and LHC energies. 
Consequently, we may neglect baryon number and as- 
sume longitudinal boost invariance ll| . We also need the 
constitutive relations for the dissipative currents. Here, 
we only consider the shear stress tensor tt^'^, the evolu- 
tion of which we describe in the approach of Israel and 
Stewart jH, {Dnf"') = ^ {2r]a'"' - ^^^'•') - ^^^^"'^xu^, 
where D — u^^d^, a^"^ = V^'^w''^, and the angular brack- 
ets <> denote the symmetrized and traceless projection, 
orthogonal to the fluid four- velocity u'^. We have also 
taken the coefficient of the last term in the massless limit. 
For details, see Ref. [l3| . 

We solve the conservation equations numerically by us- 
ing the SHASTA algorithm, see e.g. Ref. [13]. The relax- 
ation equations for the components of tt'^'^ are solved by 
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TO [fm] 


eo [GeV/fm^] 


r,„,, [MeV] 


200 


1.0 


24.0 


335 


2760 


0.6 


187.0 


506 


5500 


0.6 


240.0 


594 



TABLE I. Initialization parameters for different collisions. 

discretizing spatial gradients using centered second-order 
finite differences. We found that, in contrast to SHASTA, 
this method produces numerically stable solutions also 
for low-density matter at the edges of the system. 

With longitudinal boost invariance, we need to specify 
the values of the energy-momentum tensor in the trans- 
verse plane at some initial time tq. We assume that the 
initial energy density profile is proportional to the density 
of binary nucleon-nucleon collisions as calculated from 
the optical Glauber model (model eBC in Ref. [HI). The 
initial transverse velocity and tt^'' are set to zero. The 
maximum energy densities eo in central collisions (impact 
parameter b = 0) are chosen to reproduce the observed 
multiplicity in the 0-5% most central ^/SNN = 200 GeV 
Au-HAu collisions at RHIC El] and ^/sWN = 2.76 TeV 
Pb-|-Pb collisions at LHC [16|. For the ^/snn = 5.5 
TeV Pb-|-Pb collisions at LHC we use the multiplicity 
predicted by the minijet + saturation model [l7| . The 
initialization parameters are collected in Table HI 

Our equation of state (EoS) is a recent parametrization 
of lattice-QCD data and a hadron resonance gas [s95p- 
PCE of Ref. [18]], with chemical freeze-out at a temper- 
ature Tchcm = 150 MeV implemented as in Ref. (l9| . 

Hadron spectra are calculated by using the Cooper- 
Frye freeze-out description [20| with constant decoupling 
temperature Tdcc = 100 MeV, which will be shown below 
to give reasonable agreement with both the pT-spectrum 
and the elliptic flow coefficient for pious at RHIC. For 
the sake of simplicity, we include viscous corrections to 
the equilibrium distribution function /o as for Boltzmann 
particles, even though /o obeys the appropriate quantum 
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21|: 



f{x,p) = fo + 6f = .h 



2T2(e+p) 



(1) 



where p is pressure and p'^ is the hadron four-momentum. 
Two- and three-body decays of unstable hadrons are in- 
cluded as described in Ref. 22]. We include resonances 
up to mass 1.7 GeV. 

The shear viscosity to entropy density ratio is 
parametrized as follows. For the hadronic phase, it re- 
produces the results of Ref. [i^. In the QGP phase, rj/s 
follows the lattice QCD results of Ref. [24|. Then, rj/s 
has to assume a minimum value at a certain tempera- 
ture; in our case we take rj/s = 0.08 at T = 180 MeV. 
This is the same parametrization as used in Ref. (25| . In 
total we have four cases, see Fig. [1] (LH-LQ) rj/s = 0.08 
for all temperatures, (LH-HQ) r]/ s — 0.08 in the hadron 
gas, and above T = 180 MeV rj/s increases according to 
lattice QCD data, (HH-LQ) below T = 180 MeV, 7]/s 




0.10 



0.20 



0.30 0.40 0.50 0.60 
T [GeV] 

FIG. 1. (Color online) Different parametrizations oirj/s as a 
function of temperature. The (LH-LQ) line is shifted down- 
wards and the (HH-HQ) line upwards for better visibility. 

is that of a hadron gas, and above we set "q/s = 0.08, 
(HH-HQ) we use a realistic parametrization for both the 
hadron gas and the QGP. For the relaxation time we use 
a result motivated by kinetic theory = 577/ {e-\-p) [26[. 

Figure [5^ shows the pT-spectrum of positive pions 
in the 0-5 % most central ^/snn — 200 GeV Au-|-Au 
collisions at RHIC. Our calculations are compared to 
PHENIX data [15;]. All the different parametrizations 
of ry/s give similar agreement with the \ow-pT pion spec- 
tra. For pt ^ 1.0 GeV, the parametrizations (LH-HQ) 
and (HH-HQ) start to give slightly ff after spectra. While 
the effect of the QGP viscosity on the p^-slopes is small 
for our comparatively long initialization time tq = 1.0 
fm, it becomes more pronounced for smaller values of tq. 
On the other hand, the slopes of the spectra are almost 
independent of the hadronic viscosity and this conclusion 
remains true at least for tq = 0.2-1.0 fm. 

Figures [2Jd and [21: show the spectra for ^/snn — 2.76 
TeV and 5.5 TeV Pb+Pb collisions, respectively. Here we 
observe a much stronger dependence of the pr-spectra on 
the high-temperature values of rj/s, but the main reason 
for this is the earlier initialization time tq — 0.6 fm. On 
the other hand, the pT-spectra are independent of the 
hadronic viscosity also at LHC. 

In Figs. [2}1, [2^, and [2f we show the elliptic flow co- 
efficients for charged hadrons in the 20-30% central- 
ity class for ^/snn = 200 GeV Au-t-Au collisions and 
y/JJ^^N = 2.76 TeV and y^IWW = 5.5 TeV Pb+Pb colli- 
sions, respectively. In Fig.[2Jl the results from the hydro- 
dynamic simulations are compared to STAR 4-particle 
cumulant data :27^ and in Fig. [2]3 to recent data from 
the ALICE Collaboration 

We immediately see that, for RHIC, the four 
parametrizations for 77/ s produce values for the elliptic 
flow that fall into two classes. The curves are largely 
insensitive to the values of 77/s in the QGP phase and 
follow the value of the viscosity in the hadron gas: the 
parametrizations (LH-LQ) and (LH-HQ) with constant 
77/5 in the hadron gas result in larger V2{pt) than the 
parametrizations (HH-LQ) and (HH-HQ) with realistic 
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FIG. 2. (Color online) Transverse momentum spectra of positive pions in the 0-5% most central collisions and elliptic flow 
coefficients in the 20-30% centrality class at RHIC and LHC. Different curves correspond to the different parametrizations of 
the temperature dependence of rj/s. Data in panel (a) are from Ref. [l^ and in panels (d) and (e) from Refs. [2?!. [2^. 

rj/s in the hadron gas. We have confirmed the insensitiv- initial conditions (the same non-zero initial tt^'' for all 



ity to the values of rj/s in the high-temperature QGP 
phase by decoupling the system at Tdoc = 170 MeV. 
In that case, V2(j)t) is largely independent of the rj/s 
parametrization. The separation of curves occurs in the 
subsequent evolution in the hadronic phase. This shows 
that, within this model and at RHIC, viscous effects from 
the hadron gas dominate over viscous effects from the 
QGP, see also Refs. [i,[l3l- Due to the strong longitudinal 
expansion, the initial shear stress enhances the transverse 
pressure and thus the buildup of the flow anisotropy, 
but this is counteracted by the viscous suppression of 
anisotropics. Our simulations suggest that at RHIC these 
two effects cancel each other in the QGP phase. 

The main reason for the hadronic suppression of W2(pt) 
are the viscous corrections Sf to the particle distribu- 
tion function. Thus, the values of n'^'^ on the decou- 
pling boundary are significantly larger in the case with 
large hadronic rj/s. On the other hand, the azimuthal 
anisotropics of the hydrodynamic flow field are quite sim- 
ilar in all cases. This is demonstrated in Fig.|3l where we 
plot V2{pt) of pions at RHIC without Sf. All curves are 
much closer to each other, indicating that the space-time 
evolution in the hadron gas is similar in all four cases. 

We have tested that these conclusions arc unchanged 
if we use different tq = 0.2-1.0 fm, different EoSs, e.g. 
with or without chemical freeze-out, use non-equilibrium 



four cases), or shift the rj/s parametrizations up by a 
constant value, such that 7//s at T = 180 MeV is five 
times the AdS/CFT lower bound. Although V2{pt) and 
the slopes of the py-spectra change when we change the 
setup, the observed sensitivity of V2{pt) on the viscosity 
around T ^ 180 MeV and below, rather than on the high- 
temperature QGP viscosity is quite generic at RHIC. If 
we increase rj/s above T — 200 MeV by a factor of ten 
in parametrization (HH-LQ), the elliptic fiow is practi- 
cally the same as shown in Fig. [2Jl. This confirms that 
the value of rj/s in the high-temperature QGP phase has 
no effect on the final observable V2{pt) at RHIC, even 
though during the evolution the system spends approx- 
imately equal times above T ^ 200 MeV and between 
T - 170 and 200 MeV. 

Interestingly, the sensitivity of V2{pt) to the QGP vis- 
cosity increases with increasing collision energy, while 
the sensitivity to the hadronic viscosity decreases. This 
can be seen in Figs. ^ and [2f, which show V2{pt) f^or 
,JIWn = 2.76 TeV and = 5.5 TeV Pb-|-Pb colli- 

sions, respectively. 

At the highest LHC energy, the behavior of V2{pt) is 
completely opposite to that at RHIC. It is almost inde- 
pendent of the hadronic viscosity, but sensitive to the 
QGP viscosity. In contrast to the RHIC case, at LHC 
the differences in V2{pt) are mostly due to the difference 
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FIG. 3. (Color online) Same as Fig. [2J1, but for elliptic flow 
of thermal (i.e., without decays) pions with 5f = 0. 

in the transverse flow profiles (caused by the different 
QGP viscosities) and not due to the viscous corrections 
to the distribution function at freeze-out. The latter are 
much smaller than at RHIC: the magnitude of 5f is the 
difference between the curves (LH-LQ) and (HH-LQ) or 
(LH-HQ) and (HH-HQ) in YigM- We have also checked 
that V2{pt) at low-pr remains insensitive to the hadronic 
viscosity, even if we increase the hadronic 77/s in such way 
that it reaches 77/5 = 1.0 at T = 100 MeV, but keep the 
minimum oi rj/s fixed. The collisions at ^/snn = 2.76 
TeV are between these two extreme behaviors, as elliptic 
flow depends both on the hadronic and the QGP 77/s. 

There are several reasons why the effect of rj/s on the 
elliptic flow at LHC is so different from that at RHIC: 
first, the longer lifetime of the QGP phase, which results 
in a stronger dependence of the transverse flow on the 
viscous properties of the QGP. Second, once the system 
decouples, it has much larger transverse size and velocity 
gradients are smaller. Subsequently, dissipative effects 
from the hadronic stage are smaller and have less effect 
on the observed V2(pt)- 

In conclusion, we have investigated the effects of a 
temperature-dependent rj/s on the hadron spectra and 
elliptic flow coefficients at y^SNN = 200 GeV Au-t-Au col- 
lisions at RHIC and ^/snn = 2.76 TeV and ^/snn = 5.5 
TeV Pb+Pb collisions at LHC. We found that in all cases 
the slopes of the pt spectra of pions depend mainly on the 
high-temperature 77/s, and hardly at all on the hadronic 
viscosity. 

The effect of 77/s on the differential elliptic fiow V2{pt) 
is more subtle. At RHIC energies, V2{pt) is highly sensi- 
tive to the viscosity in hadronic matter and almost inde- 
pendent of the viscosity in the QGP phase. In contrast, 
at the highest LHC energy the opposite holds: elliptic 
flow is almost independent of the hadronic viscosity, but 
depends strongly on the QGP viscosity. Thus the extrac- 
tion of an 77/s-value for the QGP, except for its value at 
the expected minimum around Tc, is basically impossible 
using the elliptic fiow data at RHIC alone. On the other 
hand, a determination of the temperature dependence of 
?7/s in the QGP phase from elliptic flow data seems to 
be possible at LHC. This could allow the observation of 
a possible transition from the strongly coupled plasma 



near Tc, see e.g. Ref. [2^, to the weakly coupled QGP. 
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